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Abstract 

We examine tlie effect of decolierence and intermolecular interactions (chiral discrimination en- 
ergies) on the chiral stability and the distinguishability of initially pure versus mixed states in an 
open chiral system. Under a two-level approximation for a system, intermolecular interactions are 
introduced by a mean-field theory, and interaction between a system and an environment is mod- 
eled by a continuous measurement of a population difference between the two chiral states. The 
resultant equations are explored for various parameters, with emphasis on the combined effects 
of the initial condition of the system, the chiral discrimination energy and the decoherence. We 
focus on factors affecting the distinguishability as measured by population difference between the 
initially pure and mixed states and on the chiral stability as measured by the purity decay. 
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I. INTRODUCTION 



The existence of stable chiral molecules has been one of the most puzzling problems in 
quantum mechanics [1]. Since the Hamiltonian of chiral molecules is symmetric with respect 
to the parity operation, the eigenfunctions of the parity-invariant Hamiltonian should also 
be parity eigenfunctions in the absence of energy degeneracy [2], i.e., the left- and right- 
handed states of chiral molecules are not true stationary states. However, some enantiomers 
are observed to be stable for a remarkably long time despite their lack of parity symmetry. 
In the early days of quantum-mechanics Hund formulated the problem of the stability of 
chiral molecules in terms of quantum tunneling in a symmetrical double- well potential [1]. 
The two chiral states are represented by coherent superposition of the eigenstates. When 
only the lowest two eigenstates, |+) (even parity, ground) and |— ) (odd parity) with energies 
and E_, are relevant in a description of the system, the two chiral states, \L) and \R), 
localized around the potential minima, may be written as 

1^) = + (1) 

\R) = (2) 

Hence, if the system is initially prepared as one of the two chiral states, say, \R), it will 
oscillate between \R) and \L) through the inter- well barrier. The tunneling time is inversely 
proportional to the energy splitting between the two lowest eigenstates. An exceedingly 
long tunneling time corresponding to an extremely small energy splitting was thought to be 
responsible for the stability of chirality. However, upon a quantitative test, Hund's theory 
seems rather insufficient since the tunneling times observed for stable chiral molecules vary 
in a wide range, from sub-picosecond to the age of the universe, depending on the mass of 
the particle and the nature of the potential, and yet no quantum tunnelings are observed in 
these molecules. Interest in the paradox of chiral stability was redrawn by the discovery of 
weak neutral currents which violate parity symmetry at the level of nonrclativistic quantum 
mechanics, and thus may lead to stabilization of one of the two chiral configurations. How- 
ever, while theoretical calculations for this effect predict smaller energies of predominant 
(L)-amino-acids and (D)-sugars than their corresponding enantiomers, the energy difference 
is too small to be detected experimentally so far [3]. On the other hand, some have noted 
that a nonlinear Schrodinger equation does not require the eigenstates of the potential's 
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symmetry operators as the stationary states [4]. More recently, Vardi showed that the ef- 
fect of a difference between homochiral interactions (the intermolecular interactions between 
molecules of the same chirality) and heterochiral interactions (the intermolecular interac- 
tions between molecules of the opposite chirality), due to resultant nonlinearity of equations, 
can lead to far more stabilized chiral states than expected for isolated molecules [5] . Several 
authors also have pointed out that the interaction of molecules with the environment should 
be introduced to explain the stability of chiral molecules. The environmental effects on the 
system leads to decoherence, i.e. , coherence loss, and thus to a suppression of coherent tun- 
neling oscillations. Intermolecular collisions [6], interaction with photons [7], and interaction 
with phonons [8] have also been employed as the origin of the dephasing process. Another 
problem related to the issue of stabilization of chiral molecules is controlling the chirality by 
exploiting lasers: several proposals have been suggested for the preparation and detection of 
the superposed chiral states utilizing phase-controlled ultrashort laser pulses [9], and grow- 
ing attention has been directed towards using lasers to manipulate the molecular chirality 
[10]. Considering that effects of the decoherence and the intermolecular interactions on the 
system and thus on the control scenarios are far from well-understood despite their expected 
serious impacts, it is thus of great interest to understand the coherent tunneling dynamics of 
chiral systems coupled to external environments more thoroughly. In the related work [11], 
we considered an open chiral system, where a two-level approximation was taken for the 
system, and interaction between a system and an environment was modeled by a continuous 
position measurement, an approach often used in the studies of the quantum Zeno effect 
[12] and in the studies of decoherence effects in chemical reactions [13]. In particular, we 
focused on the interplay of the initial coherences and the decoherence on the distinguisha- 
bility of initially pure versus initially mixed states as measured by a population difference 
between the two chiral states, and on the vulnerability to decoherence as measured by a 
purity decay. The results provided some answers to three fundamental questions arising in 
open chiral systems as follows. 1) Can we distinguish a pure state from a mixture with same 
populations as those of the pure state? If the two chiral states are initially taken to be real 
for the pure state, no distinguishability was shown to exist for all time and any dephasing 
rate. Otherwise, the distinguishability exists, although it disappears over time faster with 
increasing the dephasing rate. 2) Why are some chiral systems stable for a very long time? 
3) Why are superpositions of stable chiral states not observed experimentally? Our result 
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suggested that the system, if it is initially a well-localized state, would tend to preserve its 
chirality when coupled to an environment that destroys left-right coherences, provided that 
the interaction time scale (dephasing time) is faster than the tunneling time scale. On the 
other hand, the system, if initially a superposition of two chiral states, would become very 
quickly a racemic mixture in the presence of the very same environment. In this paper, 
we extend the open chiral system considered in Ref. [11] to include chiral intermolecular 
interactions which are introduced within a mean-field approximation as used in Ref. [5]. 
To better understand the combined effects of the two interactions, namely the decoherence 
effect and chiral intermolecular interactions, we first consider the system in the presence of 
each interaction separately, review some results of Ref. [5], [11], and then further obtain the 
following results. 1) For a system with decoherence only and with an initially localized state, 
decreasing the tunneling frequency tends to suppress the decoherence process. It suggests 
that the more classical system is less vulnerable to decoherence. 2) For a system with chiral 
intermolecular interactions only, the nonzero chiral discrimination energies can bring about 
the distinguishability between the initially pure and the mixed states even when the chiral 
states are taken to be real for the initially pure state, where no distinguishability exists in 
the absence of chiral intermolecular interactions. It makes it clear that Vardi's claim in the 
Ref. [5] ( "• • • In the former case, spontaneous symmetry breaking may take place, turning 
a nearly racemic mixture into an optically active ensemble' ), is unjustifiable since he used 
initially pure states without any decoherence source, and thus initially pure states should 
remain pure, and initially pure state and mixed state are shown to behave differently in 
the presence of intermolecular interactions. Finally, we introduce both the decoherence and 
chiral intermolecular interactions into our model system. The dephasing time is varied to be 
very long compared with tunneling time. In this tunneling dominant region, the decoherence 
process in the absence of chiral intermolecular interactions is very slow, but, however, when 
combined with the chiral intermolecular interactions, the decoherence process tends to be 
accelerated. For a given chiral discrimination energy, we observe that, similar to the case 
with decoherence only, as the dephasing rate increases, the decoherence process becomes 
faster. On the other hand, for a given dephasing rate, the decoherence process is seen to be 
influenced by the initial condition of the system and the chiral discrimination energy. That 
is, the vulnerability of the system to decoherence depends on whether the initial state is 
strongly delocalized or localized, and whether chiral discrimination energy favors homochiral 
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mixtures or heterochiral mixtures. The organization of this work is as follows. Section II 
presents formulations of a chiral system in the presence of intermolecular interactions and 
decoherence effects. Section III gives brief descriptions of the initial conditions, measures of 
the distinguishability and the coherence loss, and computational methods. In Sec. IV the 
system in the presence of decoherence only is examined with a particular focus on how the 
initial coherence, dephasing rate, and tunneling frequency affect the distinguishability and 
the purity decay. In Sec. V the system in the presence of intermolecular interactions only is 
examined with an emphasis on the distinguishiability of the initially pure and mixed states, 
as afforded by the intermolecular interactions. In Sec. VI both effects, the decoherence and 
the intermolecular interactions are introduced, and the combined effects arc examined on 
the purity decay rate and the stationary values of the population and the coherence, with 
several different initial states for the system. Finally, Sec. VII provides conclusions and 
discuss further studies. 



II. FORMULATION: CHIRAL SYSTEMS IN THE PRESENCE OF DEPHASING 
AND INTERMOLECULAR INTERACTIONS 

To study the behavior of chiral molecules we adopt the framework of Hund's double-well 
model [1]: the chiral molecule is characterized by a symmetric, one-dimensional double-well 
potential V{x) and the x coordinate corresponds to the chiral configuration. Also we take a 
two-level approximation for the system, i.e. , consider only the two lowest energy eigenstates, 
this is valid when the energy is sufficiently low that other higher energy eigenstates are not 
involved in the system dynamics. The total wavefunction of the system may be written as 

Mt))^aL{t)\L) + aR{t)\R), (3) 

The Hamiltonian for an isolated system in the two-level approximation takes the form 

^0 = E+\+){+\ + E_\-){-l (4) 
= EM{R\ + \L){L\) - S{\L){R\ + \R){L\) (5) 

where Ej^ — ^^+~^^-^ and S — (^-~-^+) = /j^^ where u is the tunneling frequency. Substitu- 
tion of Eq.(5) into the time-dependent Schrodinger equation ih^\^) = ^o|*) leads to a set 
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of two linear equations for the expansion coefficients and ai- 

ifklL = EmttR - SttL, (6) 

ihciR = E^ttL - Sqr. (7) 

The time evolution of and shows a tunneling motion of the system between the two 
wells with oscillation frequency ou. Note that the tunneling time (= 27r/a; = 2Trh/S) is 
inversely proportional to the energy difference between the states |+) and |— ). To take 
into account the effect of the intermolecular interactions, in particular, the homochiral in- 
teractions and the heterochiral interactions, we choose Vardi's approach [5] . In Ref . [5] he 
employs a Hatrcc-Fock technique, where each molecule sees only the mean field that is gen- 
erated by the interactions of all the other molecules. The homochiral and the heterochiral 
interactions lead to two components of the self-consistent field, Vhom and Vh^t, respectively. 
Then, Eqs. (6), (7) become 

o'l = {EmdL - Sgr) - - {Vhom\(lL\^ + Vhet\aR\^) dL, (8) 

a'R = {EmttR - 5aL) - - {Vhom\aR\'^ + Vhet\(lL\'^) dR. (9) 

Here, the time-evolution is governed by two terms. The first comes from the dynamics 
of the isolated system which generates a coherent tunnehng between two wells. The sec- 
ond term involving Vhom a-nd Vhet is responsible for the intermolecular (hetero/homochiral) 
interactions between the single molecule with the other surrounding molecules within the 
mean-field approximation [5]. Note that introducing the intermolecular interactions results 
in a nonlinear Schrodinger equation of which the stationary states does not need to be eigen- 
states of the parity operator [4]. Vardi showed that if Vhom and Vhet are different, the second 
term involving nonlinear terms may affect the equation of motion significantly, and can lead 
to more stabihzed chiral states or achiral states than expected for the corresponding isolated 
molecule, depending on the initial condition and the difference between Vhom and Vhet [5]. 
However, in the mean field approximation fluctuation terms arc ignored. These fluctuation 
terms arise from the interaction between molecules and are expected to result in decoher- 
ence. Thus, in Ref. [5] the system does not go through coherence loss and never becomes 
a mixed state when starting from a pure state. Below we combine the intermolecular inter- 
actions and the decoherence effect to describe the dynamics of chiral molecules. In order to 
simulate the decoherence effect, we consider the environment surrounding the system: the 
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environment can be thought of as other degrees of freedom of the system or an inert gas or 
a condensed phase of no chirahty. We model the effects of interactions between the system 
and the environment as the continuous measurement of the system coordinate x [14, 15]. 
Then the system density operator p, which is obtained by tracing the total density operator 
over the environment, obeys the master equation [13], 

f = -^[^,P]-7'[x,[£,p]], (10) 

with 7' being a coupling strength between the system and the environment. The first term 
represents the system dynamics in the absence of the environment, and the second term, a 
multiplication of 7' and the double commutator, is responsible for the environmental effect 
on the system and is expected to destroy coherences between the eigenstates of the position 
measurement operator, x. This loss of coherence, i.e. , decoherence, has been argued as 
the essential ingredient to reconcile the two conflicting descriptions of chiral stability, the 
classical "locahzed" one and the quantum "delocalized" one [12, 16]. To invoke a two- 
level approximation for the system, we assumed that the potential barrier is very high, or 
equivalently the energy of the system is sufficiently low that the higher energy levels arc not 
involved in the dynamics. In addition, since we introduced the environment, we should point 
out that the position measurement of the system is known to increase both its energy (H) 
and energy width SE = ^ {H^) — {H^, where (•) denotes an ensemble- average, over time 
[13, 17]. The system energy and energy width will eventually become greater than the next 
higher energy levels, leading to the breakdown in the validity of the two-level approximation. 
The stronger 7', the shorter the time during which the two-level approximation will be valid. 
Using the two- level approximation and the evenness of |-|-) and oddness of |— ), we derive 
the position measurement operator as 

X w + H.C 

where H.c denotes the Hermitian conjugate and dz = {\R){R\ — \L){L\). Here ixmin are 
the positions of the double-well minima, which for a deep well are very close to the average 
positions of the right and left wave functions (denoted by ±Xavg with deflnitions Xavg = 
{R\x\R) and — Xavg = {L\x\L)). Then the master equation for the two-level system becomes 
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[18] 

f = ^[^,P]-7[^.,['^.,P]], (12) 

with the dcphasing rate 7 = 'j'x'^^^. This equation gives the evolution of the two-level 
system subjected to representing a continual measurement of the system coordinate x. 
In the first term we will use the Hamiltonian of the two-level system with the intermolecular 
interactions under a mean- field approximation as used in Eqs. (8) and (9). The second term 
in this equation represents a measurement of a population difference between two states \L) 
and \R), and leads to decay of the coherences between the \L) and \R) states (decay of the 
off-diagonal matrix elements, Plr)- To amalgamate Eqs. (8), (9) and Eq. (12), we rewrite 
Eq. (12) into the stochastic Schrodinger equation for the state vector |^) [18, 19]: 

m = (^-'-H - 7 - (^.)i*>]') mdt + V7 - (^.)i*>] mdi. (13) 

where {d'z)\^i) = (\E'|(J2|\E') and ^ is for a complex Wiener process [19]. Combining the 
intermolecular interactions and the decohcrcncc leads to the resultant Schrodinger equation 
in terms of ai, and qr, which is stochastic and nonlinear: 

ctL = --^ {EmaL - Sgr) - ^ {Vhom\(iL\'^ + VhetlaRl"^) gl - (^7l«i?l^ + \/7CWI«i?r) «l(14) 
a'R = {EmttR - Sul) - ^ {VhomlciRl'^ + VketlaL^) aR - ^7|aL|^ - VtCWIolI^) 0^(15) 

Here, the time-evolution is governed by three components. The first and the second com- 
ponents are responsible for the dynamics of the isolated system and for the intermolecular 
(hetero/homochiral) interactions between the single molecule with the other surrounding 
molecules within the mean-field approximation [5], respectively. The last component mod- 
els the interaction with the environment that induces the loss of quantum coherence in the 
system. Noting that |aRp -|- |aLp = 1, and rescaling the time as ut ^ t, Eqs. (14) and (15) 
become 

ql = -i(fl + v\aR\'^) ql + ittR - (r\aR\^ + Vff){t)\aR\'^^ ql, (16) 
a'R = -i{fl + v\aL\'^)aR + iaL- (r\aL\^ -Vrrj{t)\aL\'^^aR, (17) 

with Q = {Em + Vhom)/^, V = {Vhet " Vhom)/^, T = -f/u and drj = ^/cod^. We gradually 
build up a more complicated system starting from a simple system. In section IV we study 
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the system in the presence of decoherence effect only, using Eq. (12) with H — Hq where 
Hq is given in Eq.(5). In section V we study the system in the presence of intermolecular 
interactions only, using Eqs.(8) and (9). In section VI, we study the system in the presence of 
the combined effects of the decoherence and the intermolecular interactions, using Eqs.(16) 
and (17). 

III. THE MODEL 

A. Initial condition: Localized and delocalized states 

We examine three cases for the pure state with an initially different coherence, or equiv- 
alently with an initially different degree of localization. For an initially localized state (de- 
noted as LS), we choose \R), located in the right-well. For an initially delocalized state, we 
examine two cases, a weakly delocalized state (denoted as WDS) and a strongly delocalized 
state (denoted as SDS), with ^/oM\L) + Vol5\R) and Vo^\L) + VoIl\R), respectively 
For the LS case, the initial coherence (= plr) is 0, and the initial population difference 
between the two chiral states, Z (= pjm — pll) is 1. For the WDS case, initially Plr— 
V0.05 X 0.95 0.2, and Z=0.9. For the SDS case, initially plr^ VO.49 x 0.51 0.5, and 
Z=0.02. 

B. Measures for the distinguishability 

With the respect to the question of whether one can distinguish a pure state from a corre- 
sponding mixture with the same populations in the right and left wells. Harris and coworkers 
[21] showed, in the framework of the wavefunction, that no parity sensitive experiment can 
measure the difference between them if the states \L) and \R) are taken to be real. They 
assumed that instantaneous measurements occur on a time scale short compared to that for 
the system dynamics (tunneling) and compared to that for any interactions with the sur- 
rounding molecules. On the contrary, we include an external environment that destroys the 
right-left coherences in the framework of the density-matrix, and the intermolecular inter- 
actions (the homochiral and heterochiral interactions) in the framework of the wavefunction 
within a mean-field theory. 
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We investigate the distinguishability between the initially pure and mixed states by mea- 
suring a population in the right- well. For the measure, we compare prr of the initially pure 
vs. mixed state, and for the system with decoherence only, we use AZ(= — Z^) where 
superscripts P and M stand for the initially pure state and the corresponding mixed state. 

To examine the distinguishability between the initially pure state and the corresponding 
mixed state over time, we are required to select the initial conditions for the pure state 
{X^,Y^,Z^) and the mixed state (X^,Y^,Z^) where X = 2Re(pLR), Y = 2Im(pLR), 
and {X, Y, Z)—{Xo, Yq, Zq) ai, t — 0. For the maximum contrast between the pure and the 
mixed state, we choose the same initial population of \L) and \R) in the initial pure and 
mixed states, i.e. Zq = Z^, but choose the degree of initial coherence to be zero for the 
mixed state (maximum randomness), i.e. Xq^ = Y^^ = 0, while varying initial coherences 
for the pure state. In other words, when the initially pure state is given by 

\-^{t = 0)) = aLo\L)+aRo\R), (18) 

or, equivalent ly, 

pP(t = 0) = \aLo\''\L){L\ + \aRo\'\R){R\ + aLoa*Ro\L){R\+aRoalQ\R){L\, (19) 

the corresponding mixed state is given by 

p^{t = 0) = \aLo\^\L){L\ + \aRo\^\R){R\. (20) 

This initially mixed state may be prepared by making the system interact with an environ- 
ment that extinguishes the left-right coherence of the system sufficiently fast so that the 
system dynamics is ignorable. 

C. Measures for the degree of decoherence: Purity and off-diagonal elements 

In order to measure coherence decay, i.e. , decoherence, we use both the purity and the 
off-diagonal elements of the reduced density matrix. The purity <; is defined as [22] 

<^^Tr{p'), (21) 

where Tr denotes a trace over the system of interest. Note that the decoherence process 
does not preserve the purity of the state, that is, the purity of the resultant mixed state 
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becomes less than 1, while that of the pure state is 1. In particular, for a two- level system 
purity becomes, in the \L), \R) basis, 

Tr(p^)^plL + '2\pLR\^ + P%R, (22) 

and is given by, in the Bloch-vector representation, 

q=hx^ + Y^ + Z^ + l). (23) 

As can be easily seen in Eq. (23), 1/2 < < 1. When the purity decays to a minimum 
value, 1/2, this corresponds to X = Y = Z = 0, which represents a fully mixed state (a 
racemic mixture), with equal populations in the right and left wells and with no remaining 
coherence between them. On the other hand, when the purity remains close to 1, it implies 
that the system is resistent to the decoherence process. Another more detailed, albeit 
representation-dependent, measure of decoherence is provided by the decay of off-diagonal 
reduced density matrix elements. The diagonal elements are the probabilities of finding 
the system in the corresponding basis state. The off-diagonal elements are a measure of 
coherence. Specifically, in the \L), \R) basis, pll, Prr are the probabilities of finding the 
system in the state \L) (left well) and \R) (right well), respectively, and plr is the coherence 
between the two states. When numerically examining the decoherence rate in terms of the 
purity decay or the decay of off-diagonal elements, the initially pure states are chosen, rather 
than the initially mixed states, in order to maximize the influence of initial coherences. That 
is, <j(t = 0)=l iXi + Yi + Zi = l). 

D. Computational methods 

For the system with decoeherence only, we obtained an analytical solution [11]. However, 
for the system with intermolecular interactions only and the system with both the intermolec- 
ular interactions and the decoherence effect, we had to employ computational methods. The 
terms involved with the stochastic terms rj in Eqs.(16) and (17) are integrated by a first order 
scheme, and the other terms are integrated by the fourth-order Runge-Kutta method. The 
values of and qr resulting from the stochastic Schrodinger equation (16) and (17), when 
averaged over many realizations of the Wiener process, provide the density matrix elements, 
probabilities and coherences, according to the following relations: {L\p\L) — pll — {aLaL,): 
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{R\p\R) — prr — (a^alj), and {L\p\R) — plr — {aLa*j^. A time step of 10 ^ to 10 ^ and an 
ensemble of 10^ realizations were used. 



IV. SYSTEM WITH DECOHERENCE ONLY 



To provide a simple physical picture focused on the decoherence effect in affecting a 
tunneling motion between the two wells, we consider a system with decoherence only, i.e., 
there are no intermolecular interactions and the chiral discrimination energies are zero in 
Eqs. (16) and (17) {v — 0). In our work [11], we study the roles of initial coherences of 
the system and dephasing rates on the distinguishability of initially pure versus initially 
mixed states, and the vulnerability to decoherence, i.e. , how fast the system will lose its 
coherence or become a randomly mixed state. In this section we will brieffy review the main 
results of Ref. [11] and then further discuss the effect of the tunneling frequency u) on the 
distinguishability and the purity decay rate. 



A. Review: Effect of initial coherence and dephasing rate 

Within a two-level approximation, using H — Hq, for simplicity choosing the energy 
origin to correspond to Em = 0, i.e. , = —E_, and switching from the density- matrix 

elements of the states \L) and |-R), pll, Prr, and plr to the real quantities, X, Y, and Z 
for the Bloch-vector representation, Eq. (12) gives the following set of equations (see, for 
details, Ref. [11]) : 

^ = -47^> (24) 
dY 

_ ^ -4jY + 2ujZ, (25) 

f = -^^y- (26) 

The resultant equations are of standard form [23]. At long times, when the system reaches 
equilibrium, dX/dt = dY/dt = dZ/d.t = 0, and the steady-state solution to Eqs. (24)-(26) 
becomes X — Y — Z — 0. This stationary state represents a fully mixed state, with equal 
populations in the right and left wells and with no remaining coherence between them, and 
is reached regardless of the initial conditions, the system's tunneling frequency ou, and the 
dephasing rate 7 (related to the nature of the environment and the potential of the system) . 
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With an arbitrary initial state [X — Xo,Y — Yq, Z — Zq], the general analytical solution (not 
shown here) can be obtained separately for 7 < a; (tunneling dominant region) and for 7 > a; 
(dephasing dominant region). The general solution reveals that, in the dephasing dominant 
regime for an initially localized state, say \R) (note that this corresponds to an initially chiral 
state), the system approaches the steady state (fully mixed state) monotonically without 
any oscillations present but, however, as 7 increases, tends to be stuck in the right well 
on a time scale of I/7 (this behavior is equivalent to the quantum Zeno effect [18]). Ref. 
[11] demonstrated that, as 7 increases in this dephasing dominant region, an initially chiral 
system tends to preserve the chirality, while an initially very delocalized state such as the 
SDS case tends to become a raccmic mixture more quickly. From a general solution to 
Eqs.(24)-(26), the difference between the time-evolution of Z of the initially pure state and 
that of the corresponding mixed state, AZ(= Z^ — Z^) is given by [11]: 



is crucial to the distinguishability. If — 0, AZ{t) — 0, which means no distinguishabihty, 
for all time and for all 7. However, if Y^^ 7^ the distinguishabihty does exist, although it 
disappears over time because the decoherence impels both the initially pure state and the 
corresponding mixture to the maximally randomized stationary state. For the tunneling 
dominant region, as 7 increases, the period of the oscillation increases and the amplitude 
decays faster, leading to AZ approaching zero more rapidly, at which point one can not tell 
whether the system was initially in the pure state or in the corresponding mixture upon the 
measurement of population. However, for the dephasing dominant region, as 7 increases, 
the distinguishability disappears more slowly, suggesting that considerably strong dephasing 
may unexpectedly afford means to observe the long-time surviving distinguishability. Noting 
that for a two-level system purity is given by Eq.(23) and using Eqs.(24)-(26), one can obtain 



In our previous work [11], the interplay of the initial coherence and the dephasing rate on the 
purity decay, in particular, as it relates to the racemization and chirality stabilization, was 
clarified. For very short times, increasing the dephasing rate 7 and the initial coherences 




with s — 2^ up' — 7^ and s — 2^ 7^ — uj"^. At a glance, one can easily see that the value of 1^ 



dt 



47(^2 + y^). 



(29) 
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{Xq + Yq — 4|piRp) increases the purity decay (decoherence) rate. On the other hand, 
for later times, for larger initial coherences {Xq + Yq 1), the purity still tends to decay 
faster with increasing 7, but for the smaller initial coherences (Xq + Yq ^ 0), the purity 
shows the opposite behavior when 7 is increased beyond a certain threshold, i.e. , the purity 
decay starts to be suppressed. This interesting observation enables us to suggest that, in 
the environment destroying the left-right coherences, on a timescale faster than the intra- 
system dynamics (the tunneling time scale) , the initially strongly localized state would tend 
to preserve the chirality, while, on the other hand, an initially strongly delocalized state such 
as a superposition of the two chiral states, would become a racemic mixture very quickly. 
It is notable that the observation agrees with the prediction that if the system dynamics 
is ignorablc compared with the interaction between the system and the environment, the 
coherence between the eigenstates of the decohering operator will be lost fastest [12, 16]. 
The applicability of the prediction can be seen by realizing that the two chiral states \L) and 
\R) correspond to the eigenstates of the decohering operator {az — \R){R\ — |-^)(-^|), and the 
stronger dephasing rate beyond a certain threshold allows us to ignore the system dynamics 
compared with the decoherence process. Then, according to the prediction, a superposition 
of the two chiral states will lose its coherence much faster than a localized chiral state (|L) or 
\R)). However, it is also notable that there was previously no prediction for the observation 
in Ref. [11] that increasing the initial coherences {Xq + Yq) leads to a system that exhibits 
more rapid decoherence. 

B. Tunneling frequency uo {=6/h) 

The typical behavior of AZ{t) for several values of u for a given value of 7 (here chosen 
as 1) is shown in Fig. 1 (a) and (b). For a; < 7 (dephasing dominant region, here u < 
1, shown in (a)), as time increases AZ/Yq first decreases from zero to a minimum and 
then increases monotonically towards zero. Note that, as uj increases, the maximum of 
\AZ/Yq\ gets larger and is obtained at an earlier time, i.e. , imax gets smaller and the 
distinguishability disappears more rapidly beyond tmax (i-©- ^Z/Yq goes to zero more 
quickly). For instance, compare the cases cu = 0.1 and u = 0.8 in Fig. 1 (a). For the 
case UJ = 0.1, IAZ/Fq^I = 0.0494 at t — imax ~ 1-5, and, is still nonzero at t = 8 having a 
value of 0.0464 which is roughly the same as the maximum value of \AZ/Yq\), indicating 
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the long-time surviving distinguishability, while for the case uj — 0.8, |AZ/yo^| = 0.315 at 
i — ^max = 0.58, and is almost zero at t = 8, indicating no distinguishability for i > 8. 
This observation suggests that a very slow tunneling oscillation may unexpectedly allow 
us to observe the long-time surviving distinguishability. On the other hand, for > 7 
(tunneling dominant region, here > 1, shown in (b)), as uj increases, the period of the 
oscillation in AZ decreases, but the amplitude envelope decays exponentially with a rate 
that is almost independent of io. In other words, for the tunneling dominant region, the 
rate with which AZ approaches toward zero is largely determined by 7, as expected from 
Eq.(27). For all the cases a; — 2,4,6 [Fig. 1 (b)], upon the measurement of population 
after t = 3, one can not tell whether the system was initially in the pure state or in the 
corresponding mixture. Figure 1(c) captures well the effect of increasing the tunneling 
frequency of the system on the purity decay at a given 7 (= 1). Initially the system was 
prepared in the \K) state (the LS case). Interestingly and non- intuitively, increasing the 
tunneling frequency uj tends to increase the rate of purity decay rather than reducing it. For 
the case a;=0.2, 0.4 and 0.8 (the dephasing dominant region), increasing a; clearly accelerates 
the purity decay for all time. Even for the case cu = 2, 4, and 8 (the tunneling dominant 
region) where the oscillatory decay becomes visible, the purity decay rate increases with 
increasing u at short times, while the overall decay timescale appears to be dominated 
by the dephasing rate 7. In regard to the quantum-classical transition, realizing that a 
smaller lo corresponds to a higher inter-well barrier, a longer tunneling time, and thus a 
more classical system, the observation that decreasing lo decreases the decoherence rate can 
be put in the following way; a more classical (quantum) system is less (more) vulnerable to 
decoherence. It is also noteworthy that the master equation of the system of interest in this 
section (the chiral system within a two-level approximation in the presence of an environment 
destroying the coherence between the two states, see Eqs. (24) to (26)) is identical to that 
resulting from a two-level atom driven by a field undergoing a pure dephasing [24] in quantum 
optics, as previously mentioned in Ref. [18]. This equivalency emerges where the following 
replacements are made: \L,K) the two energy eigenstates, io Rabi- frequency, and 
7 — 7- pure dephasing rate. Noting this equivalency enables us to recognize that the above 
observation (increasing a; tends to increase the purity decay rate) is also equivalent to the 
result of Ref. [25], where aspects of coherence and decoherence were studied with a two- 
level system interacting with a decohering environment and a resonant continuous-wave 
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electromagnetic field, that increasing the Rabi frequency increases the decoherence rate. 

V. SYSTEM WITH INTERMOLECULAR INTERACTIONS ONLY 

The system with intermolecular interactions within a mean-field approximation was stud- 
ied in Ref. [5], with a particular focus on the role of intermolecular interactions in estab- 
lishing chiral stability. Localized chiral states were shown to be stable when the homochiral 
interactions are energetically favorable to heterochiral interactions {v > 0), whereas delo- 
calized states are dominant when the heterochiral interactions are energetically favorable 
{v < 0). For a self-completeness and a later comparison with the case in the presence of 
both the intermolecular interactions and the decoherence effect, we reproduce some of the 
results of Ref. [5], with emphasis on the interplay of the chiral discrimination energy param- 
eter V and the initial condition of the system in establishing the anharmonic oscillations in 
population of the two wells of the system. In subsection V B we show that a nonzero chiral 
discrimination energy can bring about the distinguishability between the initially pure and 
the corresponding mixed states even for the Yq — case, where there is no distinguishability 
in the system in the absence of the intermolecular interactions {v = 0), regardless of the 
decoherence effect, as shown subsection IV. 

A. Review 

The system with intermolecular interactions only is described by Eqs.(16) and (17) with 
r = 0. Below for all the numerical calculations Q is chosen as 1. The time-evolution of the 
initially localized state \R) (the LS case) with a variation of v is shown in Fig. 2. For the v — 
case we see coherent tunneling from one well to the other. As |v| increases (the chiral energy 
difference \ Vhet — yhom\ increases), the nonlinearity of the differential equations increases and 
the oscillations in population between the two wells becomes anharmonic. For a sufficiently 
large \v\{> 4) and with increasing \v\, population-trapping becomes stronger and maximum 
coherence decreases. We observed a similar relation between the coherence \plr\'^ and the 
self-trapping in the WDS and the SDS cases (not shown here): when self-trapping due 
to intermolecular interactions is manifested in pnn, \plr\'^ is seen to exhibit a decrease in 
the oscillation amplitude, and thus a decrease in the maximum coherence, which prohibits 
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a complete tunneling between the two wells, and consequently results in the self-trapping 
of population. Considering that characteristic values of the chiral discrimination energy 
difference, \Vhet — Vhom\, ^-re far greater than the tunneling splitting S, i.e. , |v| ^ 1 [26], Fig. 
2(a) demonstrates that much of the system in a condensed phase can be self-trapped in one 
of the enantiometric forms to a degree which depends on the chiral discrimination energies 
(i.e. Vhet and Vhom)- As demonstrated in Fig. 2, behaviors of prr, |pLi?P, and thus the self- 
trapping effect are independent of the sign of li^l for the initially localized state (LS case). 
Note that the LS case has no initial coherence between \L) and \R), i.e. , pLR{t = 0) = 0. 
However, when there exist nonzero initial coherences, the stability of the initial states does 
depend on the sign of |f |, i.e. , whether the homochiral mixture is energetically favorable 
{v > 0) or the heterochiral mixture is favorable {v < 0). We present two cases of different 
initial coherences, the WDS and SDS cases, for several values of v. For the WDS case with 
positive V [Fig. 3(a)], as v increases from to 2 the oscillation period of prr increases and 
then decreases beyond v — 2 with the population becoming self-trapped ior v > 5. This 
observation is similar to that of the LS case. However, when v is negative [Fig. 3(b)], this 
initially small degree of delocalization leads the system to oscillate between the two wells 
(chiral activity) with a faster oscillation frequency than that of the v = case, until v 
reaches quite a large magnitude, for instance, v — —10, when the system is now self-trapped 
again. Comparisons of the v — 5 and the v — —5 cases show how the initially weakly 
delocahzed state (the WDS case) can behave very differently depending on the sign of v. 
For large negative v (e.g., v — —10) the system becomes self-trapped. For large positive v 
(e.g., V = +5) the system also becomes self-trapped, but such a trapping is stable in sense 
that it remains self-trapped as v is further increased. For positive v in the WDS case the 
self-trapping becomes very stable, while for the negative v the tunneling is more assisted 
and stabilized by the addition of the intermolecular interactions. This behavior is shown 
in the exactly opposite way for the SDS case, since the initial conditions with almost equal 
populations in each well (ol ~ qr) correspond to stable states for negative v but unstable 
states for positive v. With positive v [Fig. 4(a)] the initially strongly delocahzed state 
exhibits oscillations about prr ^ 0.75 for f = 5 and prr ^ 0.7 for v = 10, indicating that 
the system is chiral active, whereas as for v = and v = 2, prr exhibits a small amplitude 
oscillations about ~ 0.5, indicating that the system is not chilar active. On the other hand, 
for negative v [Fig. 4(b)] the delocahzed state becomes very stable around prr — 0.5 with 
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an oscillation frequency that increases with increasing \v\. 
B. Distinguishability of pure and mixed states 

One of the main motivations of this paper was to answer some questions arising from 
Vardi's claim in the last paragraph of Ref. [5]. He argued that from the results "■ ■ ■ In 
the former case, spontaneous symmetry breaking may take place, turning a nearly racemic 
mixture into an optically active ensemble." Since his model does not associate with any 
decoherence source, the system remains pure if it is initially pure. All of his results were 
obtained for initially pure states. Thus when he used "racemic mixture" , it would have been 
correctly described as a delocalized pure state as opposed to a racemic mixture that is found 
in nature and associated with loss of coherences. At this point we face two questions. First, 
if the system starts as a mixed state, then will it show results similar to those of initially 
pure states? This question may be put as follows: can we distinguish the initially pure 
and the initially mixed states for the system with intermolecular interactions? If we cannot 
do so by measuring populations, for instance, prr, then Vardi's claim may be justifiable. 
Otherwise, decoherence should be seriously considered to answer the second question: how 
can an initially pure state end up as a nearly racemic mixture or an optically active ensemble 
in real system? Wc attempt to answer the first question here and the second one in the next 
section. We start by showing the time-evolutions of population and coherence in Fig. 5 with 
the initially pure and the corresponding mixed state for the WDS and SDS cases (note that 
for both of these cases — 0) in the absence of the intermolecular interactions {v — 0). 
Figure 5(a) demonstrates that when Yq = 0, there is no distinguishability in terms of 
the population {AZ = 0), as already seen in Eqs. (27) and (28), whereas there arc clear 
differences in terms of the coherence as shown in Fig. 5(b). Now we proceed the issue of 
distinguishability of the initially pure and mixed states for the system with the intermolecular 
interactions. In the previous section VA, changing the initial condition from the LS case 
to the WDS and SDS cases, and thus creating nonzero initial coherence is seen to lead the 
system to behave differently depending on the sign of v. Thus, one may imagine that the 
initially mixed state without initial coherences may behave differently from the initially pure 
state with nonzero initial coherence, even when Yq^ = such as in the WDS and the SDS 
cases. We compare the initially pure state for the WDS and SDS with their corresponding 
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mixed states, in terms of the population and the coherence, in particular, ior v — 5 and 
V — —5. Figure 5(c) demonstrates that for the WDS case a nonzero v enables one to tell 
whether the system was initially in a pure or a mixed state by measuring populations. The 
pure state is observed to behave differently, depending on the sign of v, while the behavior 
of the corresponding mixed state is independent of the sign of v due to the absence of its 
initial coherence. Since the WDS case is slightly delocalized from the LS case, the behavior 
of the population (as measured by prr) of the initially mixed state for the WDS case is 
quite close to that of the pure state of the LS case [Fig. 5(c)]. For the SDS case, the similar 
distinguishability, brought about by the nonzero v, is seen in Fig. 5(d). Let us clarify the 
behavior of the initially mixed state related to the LS case with the same v{= ±5): one can 
note that of the mixed state for the WDS and SDS cases is seen to oscillate in the same 
way as that of the initially pure state for the LS case with the same v and differences are 
found only in the amplitude and the associated mean value prr (as obtained by averaging 
over an oscillation cycle) [Figs. 5(c) and 5(d)]. When the initially pure state is given by Eq. 
(18), the corresponding mixed state is chosen by Eq. (20). Then, the probability of being in 
the right well for the initially mixed state, p^^^(t) = {R\p^{t)\R) can be expressed in terms 
of the probability of being in the right well for the LS case (denoted as |a^^(t)P): 

PRRit) = {\an\'), (30) 
= \aLoWR{t)\' + \anoWR''{t)\', (31) 
= \aLo\' + {\aRo\'-\aLo\')\af{t)\'. (32) 

where |a^^(t)P ( |fl^^(^)P) denotes the probability of being in the right well at time t for the 
state initially localized in |L) (|i?)). Here we used |a^^(t)p = 1 - |a|;^(i)p = 1 - \a^^{t)\'^: 
the second equahty holds since Eqs. (16) and (17) with F = (decoherence-free case) are L- 
R symmetric. The above Eq. (32) indicates that the probability of the initially mixed state 
should have the same oscillation period as that of the initially pure state localized in \R) (LS 
case) with the same v, but the amplitude and average of the probability of the former are 
changed from those of the latter, depending on its |aLoP; and la^joP- For the SDS case, for 
instance, since la^oP = 0.49, and la^op — 0.51, of the initially mixed state with v — ±5 
oscillates from 0.51 to 0.506 in tune with of the initially pure state locahzed in \R) (LS 
case) (= |aR^(^)P) with the same v{— ±5) which oscillates from 1 to 0.8 [Fig. 5(d)]. The 
similar observation is made in the WDS case as well [Fig. 5(c)]. As clearly demonstrated 
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in Figs. 5(c) and 5(d), for a system with intermolecular interactions, the initially mixed 
state behaves differently from the initially pure state as measured by the probabilities. In 
addition, the initially mixed state does not behave differently depending on the sign of 
V. Thus, Vardi's claim is unjustifiable. For instance, an initially nearly racemic mixture 
remains as a racemic mixture, regardless of whether the chiral discrimination energy favors 
the homochiral mixture or the heterochiral mixture, as shown in Fig. 5(d). Now it seems 
very critical to include the initial coherences adequately to obtain the correct behavior of 
the system subject to the intermolecular interactions. Further, if one were to create a nearly 
racemic mixture or an optically active ensemble in laboratory, starting from an initially 
pure state, it would be essential to include the decoherence process in the assessment and 
modeling of the system dynamics. We explore this in the following section. 



VI. SYSTEM IN THE PRESENCE OF BOTH DECOHERENCE AND INTER- 
MOLECULAR INTERACTIONS 

In this section, we finally consider the system in the presence of both a decoherence effect 
and intermolecular interactions. In Section IV we discussed general features of the decoher- 
ence effect on the system without chiral discrimination energies, i.e. , v=0. Two prominent 
system behaviours were observed: suppressing the tunnelings (oscillations) between the two 
chiral states by destroying the coherence, and forcing the system to the long-time mixed state 
with a complete randomness, p{t — >■ oo) = |(|L)(L| -|- \R){R\), regardless of the dephasing 
rate and the initial condition as long as the decoherence is present. In the study of the de- 
coherence effect on a system having the intermolecular interactions, the former is observed, 
while the latter docs not occur for all cases. The long-time stationary state is observed to 
be dependent on the decoherence strength F, the initial condition of the system, and the 
chiral discrimination energy parameter v. Details are as follows. In particular, we examine 
three different classes of initial conditions: the LS, WDS, and SDS cases. We apply Eqs.(16) 
and (17). Coherence losses are monitored by measuring the purity rather than Ipl^P since 
the purity enables us to compare easily the cases with different maximum coherences and 
is representation-independent. We point out that all decoherence-free cases were presented 
up to t = 5 (Figs. 2 to 4), while the cases with the decoherence will be presented up to 
i = 20 to reveal the longer-time behavior. We varied F from 1.6x10""^ to 1.6x10"^, which 
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means that dephasing time is very long compared with tunnehng time (tunnehng dominant 
region). Figures 6(a) and 6(b) show behaviors of the population and the purity for several 
values of f at a very small decoherence (here chosen as F = 1.6 x 10~^) for the LS case. 
None of the cases examined here reached a stationary state within t = 20. As expected, 
introducing decoherence results in decay of the oscillation amplitudes of both the population 
and the purity. Details are affected by \v\: for such a small decoherence (F = 1.6 x 10"'^), 
the V — case is seen to make no significant changes within a time window of i = 20, while 
cases with a nonzero v tend to experience relatively significant changes, suggesting that the 
intermolecular interactions make the system more vulnerable to decoherence. However, the 
decoherence rate is not exactly proportional to |f |. A relative sensitivity to decoherence 
depending on the magnitude of v is more clearly seen from comparing the purity decay 
for several cases with different v in Fig. 6(b). The order in the rate of purity decay R is 
R{v = ±4) > R{v = ±3.9) > R{v = ±3) > R{v = ±5) > R{v = ±10). From the behavior 
of prr for the decoherence-free case shown in Fig. 2(a), one can see that until li^l becomes 
large enough to cause the self-trapping (v — ±5, ±10), the purity decay rate R increases 
with increasing \v\. However, when the self-trapping occurs, R becomes smaller than other 
cases that show no self-trapping, and tends to be smaller with a stronger self-trapping, i.e. 
, greater v. This trend is also observed in the LS case with larger F, F = 1.6 x 10~^ in 
Fig. 6(c) and F = 1.6 x 10~^ in Fig. 7(a). For instance, consider the case with v — ±10, 
which induces a very strong self-trapping [Fig. 2(a)]. The purity decay rate for = ±10 is 
almost overlapped with that for v = with F = 1.6 x 10"^ [Fig. 6(b)] and F = 1.6 x 10"=^ 
[Fig. 6(c)], and even becomes smaller than that for v = with F = 1.6 x 10~^ [Fig. 7(a)]. 
It suggests that for some chiral discrimination energies, especially those inducing a strong 
self-trapping (localization) , the resistance to decoherence can be stronger than for the v — 
case, whereas, for values of v not inducing strong self-trapping the opposite effect is seen 
with respect to decoherence, namely the non-zero v case is more sensitive to decoherence 
than the v — case. It is noteworthy that the effect of chiral interactions on the LS case is 
very similar to that of decoherence; for the initially localized state, namely purity decay is 
accelerated with increasing decoherence, but starts to be suppressed with continued increase 
in decoherence beyond a certain threshold decoherence [11]. Figures 7(a) and 7(b) show that 
for F = 1.6 X 10~^ both of the purity {t oo) and pRR{t — >■ oo) for the LS case tend to 
be larger for a larger \v\, implying that for a bigger chiral discrimination energy the system 
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relaxes to a stationary state with greater chirality. Specifically, within t — 20, for v — ±3, 
the system will be relaxed to a racemic mixture {prr — > 0.5), whereas, for v — ±10, the 
system will preserve its chirahty as prr > 0.95 at i = 20 [Fig. 7(b)]. Also it is notable that 
the purity for v = ±3 has a stationary value of roughly 0.57, whereas the purity for v = ±10 
does not decay much from the initial value (=1), being 0.96 at t = 20 [Fig. 7(a)]. This can 
be understood by realizing that, from comparison with the decoherence-free cases shown in 
Fig. 2(a), V — ±10 stabilizes the initially localized state by self-trapping it to the initially 
located well, and this effect will not be much affected by the decoherence inhibiting the oscil- 
lations between the two wells. On the other hand v — ±3 does not bring about that kind of 
self-trapping effect, and rather makes the system anharmonically oscillate between the two 
wells, which would be prohibited by the loss of coherences in the presence of decoherence. 
For the LS case, the stronger the self-trapping, the less vulnerability to decoherence. Also, 
in relation to stationary values of the purity and prr, for v — ±3, the observation that 
the purity — > 0.57 and prr — > 0.5, implies that, from Eq. (22), plr does not relax to zero 
(|pLi?P goes to roughly 0.04), and the system can have therefore gained coherence due to 
the combined influence of the decoherence, the intermolecular interactions, and the initial 
condition. It is also notable that decoherence brings the purity and population {prr) for 
V — ±3.9 and ±4 cases quite close together, as shown in Fig. 7(a) and [Fig. 2(b)], respec- 
tively, smoothing out the dramatic differences between them in the absence of decoherence 
[Fig. 2(a)]. To see the effect of increasing decoherence on the population relaxation and the 
purity decay for a given v, we present prr and the purity ior v — 5, with several values of 
r, in Fig. 8. As decoherence F increases, the rate of relaxations in both the population and 
the purity increases. Also it is notable that both of PrrII — )■ oo) and the purity (t — )• oo) 
tend to decrease with increasing F, implying that with increasing F the v = 5 case tends to 
approach a less chiral stationary state. Therefore, by inducing the loss of coherences between 
the two states and thus leading to the suppression of the oscillation (tunneling) , decoherence 
can bring to a system a stabihzation of the chirality with F = 1.6 x 10~^, or a less chiral 
state (or possibly racemization) with F = 1.6 x 10~^, as determined by the interplay of the 
decoherence (F) and the intermolecular interactions (v). Now let us turn to the opposite 
case, the initially strongly dclocalizcd state (SDS) case. The results of the SDS case are 
presented in Fig. 9 with a variation of v at F = 1.6 x 10~^. First, in regard to the relative 
order of the purity decay, one can easily see that the stronger the chiral discrimination en- 
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ergy (v), the faster the purity decays, i.e. , the larger the decoherence rate is [Fig. 9]: with 
positive V, R(v = 2) < R(v = 3) < R(v = 5) < R(v = 10) [Fig. 9(b)], and with negative 
V, R{v = -2) < R{v - -3) < R{v = -5) < R{v = -10) [Fig. 9(d)]. Unhke the LS case, 
an exception is not found in the SDS case due to non-existence of self-trapping state in the 
SDS case for parameter ranges considered here [Fig. 4]. Second, the SDS case with positive 
V is seen to be more vulnerable to decoherence than that with negative v. For the SDS 
shown in Fig. 4, time-evolution of population of the initially strongly delocalized 
state (ol an) depends strongly on the sign of v. When there is no decoherence effect, for 
negative v, the oscillates in a narrow range about ~ 0.5 whereas for positive v, prr 
oscillates over a much wider range between p^R ~ 0-5 and p^i^ ^ 1. That is, the initially 
strongly delocalized state becomes stable for negative v but unstable for positive v. This 
feature continues to be effective upon introducing the decoherence effect, and thus the SDS 
case shows a faster purity decay and a faster population relaxation for positive v than for 
negative v [Fig. 9] . It suggests that a critical difference in the stabihty of the initial states, 
depending on the sign of v, i.e. , whether the heterochiral interactions are stronger or weaker 
than the homochiral interactions, affects the sensitivity of the system to the decoherence 
effects. This is reconfirmed in the WDS case (a^, ^ aji) for which results arc presented in 
Fig. 10 with a variation oi v a,t T — 1.6 x 10~^. In contrast to the SDS case [Fig. 9], the 
WDS case with negative v is more vulnerable to decoherence than that with positive v. For 
the WDS case in the absence of decoherence, when the homochiral mixture is preferred (the 
positive v), the initially weakly delocalized state becomes stabilized, while, when the hete- 
rochiral mixture is preferred (negative v), the initial state becomes highly unstable and the 
system leads to the facihtated tunnehng [Fig. 3]. Thus, upon introducing the decoherence, 
the system with positive v shows a less vulnerability to the decoherence, i.e. , a slower purity 
decay and a slower population relaxation, than the corresponding cases with negative v [Fig. 
10]. It is very interesting to see that the character of intermolecular interactions can change 
the sensitivity of the system to the decoherence, by (de)stabihzing the initial state. On the 
other hand, for the WDS case, the relative sensitivity to the decoherence seems to have a 
rather complicated dependence on v due to its medium sized initial coherence, compared 
with the two extreme cases, the LS and SDS cases. The effects of increasing decoherence 
on the WDS and SDS cases for a given v are demonstrated in Figs. 11 and 12, respectively. 
Similarly to the LS case, for a given v, in both of the WDS and SDS cases, one can see 
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that the larger the decoherence is, the faster the purity decays and the faster the population 
relaxes. Further, related to the homochiral/heterochiral preferences in the intermolecular 
interactions, introducing decoherence does not necessarily lead to the breakdown of the ho- 
mochiral/heterochiral preference. For the WDS case, for instance, a value of f = 5 yields 
a system with a homochiral preference when F = 1.6 x 10~^ {pRR{t — )■ oo) is expected to 
be greater than 0.9 as shown in Fig. 11(a)), while a value of v — —5 results in a racemic 
mixture without chirality, reflecting its heterochiral preference for the same value of F [Fig. 
11(c)]. Similarly, for the SDS case, at all F examined here, for i> = 5 prr — > 0.51 [Fig. 
12(a)], implying a small homochiral preference. For v — —5 prr is expected to relax to 0.5 
with oscillations around 0.5, implying the heterochiral preferences [Fig. 12(c)]. 

VII. CONCLUSION 

We studied the effects of chiral intermolecular interactions and decoherence on the dis- 
tinguishability between the initially pure and mixed states, and the chirality stabilization in 
a chiral system. A two-level approximation was adopted for the system. Chiral intermolec- 
ular interactions were introduced within a mean-fleld theory. The interaction of the system 
with an environment (to take account for externally induced decoherence) was modeled as 
a continuous position measurement of the system, which reduces to a measurement of a 
population difference between the two chiral states. To better understand the combined 
effects of the two interactions, we flrst investigated the system in the presence of each inter- 
action separately. For the system with decoherence only, most features is studied in Ref.[ll] 
where we obtained the following results (which we presented in summarized form in this 
work). When Yf^ is nonzero, for the dephasing dominant region, increasing the tunneling 
frequency can lead to long-time survival of distinguishability, while for the tunneling dom- 
inant region, the timescale of the distinguishability is largely dominated by the dephasing 
rate. For the initially \R) state (LS), decreasing the tunnehng frequency tends to suppress 
the purity decay, suggesting that the more classical system is less vulnerable to decoherence. 
In the system with chiral intermolecular interactions only, we showed that the nonzero chi- 
ral discrimination energies can bring about the distinguishability between the initially pure 
and the corresponding mixed states even when Yq = 0, where no distinguishability exists 
in the absence of the intermolecular interactions. Finally, we introduced both decoherence 



24 



and intermolecular interactions to the system, and then explored the combined effects on 
the purity decay rate and the stationary states of the system with several different initial 
conditions. The dephasing time was varied to be relatively very long compared with tun- 
neling time (tunneling dominant region) so that the decoherence process in the absence of 
chiral intermolecular interactions is very slow. However, when combined with the chiral 
intermolecular interactions, the decoherence process tends to be accelerated. For a given 
chiral discrimination energy v, we showed that, similar to the case with decoherence only, 
the purity decays faster as the dephasing rate increases. For a given dephasing rate, the 
purity decay was shown to be influenced by v and the initial state of the system. As the 
chiral discrimination energy |f | increases, the SDS case is observed to have a faster purity 
decay, while the LS case to have a purity decay suppressed beyond a certain threshold \v\ 
because of onset of induced self-trapping. This is very similar to the interplay of the initial 
coherence of the system and the dephasing rate on the purity decay in the system with 
decoherence only. Also we clarified the role of the sign of v associated with the initial state 
on the purity decay. For the WDS case, the purity is seen to decay slower with positive 
V than with negative v since positive v, which prefers homochiral mixtures, stabilizes the 
initially weakly dclocalizcd state. For the SDS case, the opposite is true. The results should 
serve as a prototype for understanding the results of a chiral system in more complicated 
and realistic conditions, such as systems with more than two degrees of freedom (beyond a 
two-level approximation) and with other types of environment, for instance, which induces 
both population loss (energy relaxation) and phase loss (dephasing). This work may be 
extended to an open chiral system interacting with lasers. In recent years, there has been 
increasing interest on the control of molecular chirality [9, 10]. However, the decoherence 
effect and the intermolecular interactions on these preparation/control scenarios are far from 
well understood. Realizing that the decoherence process can be detrimental to the control 
scenarios by destroying the very coherence that is needed for the control scenario, it would 
be critical to understand and include the decoherence process properly when the time scales 
of system dynamics, and/or environmental effects (energy loss, dephasing) are comparable 
to the interaction time scale between the system and the laser field. In addition, our results 
suggest that taking into account these effects may be important even when the decoherence 
time scale is not comparable to or longer than other time scales, since a system with nonzero 
chiral discrimination energies tends to be more vulnerable to decoherence, i.e. , tends to have 
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a faster purity decay than a system without chiral intermolecular interactions. 
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FIGURE CAPTIONS 



Fig. 1. Dependence of distinguishability and decoherence process on tunneling frequency 
a; as a function of time for the decoherence only case. Distinguishability versus time as 
measured by AZ{t){= Z^{t) — Z^{t))/YQ is depicted for the dephasing dominant region 
a; < 1 in panel (a) and for tunneling dominant region cu > 1 in panel (b). The values of co 
are indicated in the legend and included in both panels to provide a comparison. Purity vs. 
time is depicted in panel (c) for a range of cu values; 7 is chosen as 1. All the variables are 
in dimensionless units. 

Fig. 2. Time evolution of population and coherence |pLiiP are shown for the LS case 

with V = (thick line), ±3 (•••), ±3.9 ( ), ±4( ), ±5(- . - . -), and ±10 (thin 

hue) . 

Fig. 3. Time evolution of prr is shown for the WDS case with various values oi v. (a): 

^; = (thick line), 2 ( ), 3 (■ ■ ■), 5(- . - . -), and 10 (thin line), (b): t; = (thick line), 

-2 ( ), -3 (■ ■ ■), -5(- .-.-), and -10 (thin line). 

Fig. 4. Time evolution of prr is shown for the SDS case with various values of v. (a): f = 

(thick line), 2 ( ), 3 (• ■ ■), 5(- . - . -), and 10 (thin line), (b): v = (thick line), -10 

(thin line). Note the difference in the scale of ordinate in panels (a) and (b). 

Fig. 5. Time evolution of prr and IplrI"^ are shown with various initial conditions. For (a) 
and (h), V — and four different initial conditions are presented: the initially pure state for 

the WDS case ( ), the initially pure state (• • •) for the SDS case, and the corresponding 

mixed states are presented as (squares) and (circles), respectively. For comparison, the LS 

case with v = is also presented ( ). Panels (c) and (d) display the WDS and SDS 

cases, respectively, with four different conditions: the initially pure states with v = 5 ( ) 

and V = —5 {■ ■ ■) and the corresponding mixed states as (squares) and (xxx), respectively. 

For comparison, the LS case with v — 5, is also presented ( ). The inset in (d) is for 

the SDS case with the initially mixed state ior v — 5 (squares) with enlarged ordinate which 
reveals that it oscillates in tune with the LS case with v — 5. 
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Fig. 6. The LS case with various v with T = 1.6 x 10"^ for (a),(b) and with T = 1.6 x 10"^ 

for (c). In (a), prr vs. the rescaled time with v — (thick hne), ±3.9 (—•—), ±4( ), 

and ±5(- . - . -). In (b) and (c), purity vs. the rescaled time with v — (thick line), ±3 

(■ ■ ■), ±3.9 ( ), ±4( ), ±5(- .-.-), and ±10 (thin line). Note that the case with 

V = ±10 overlapps with the v = case. 

Fig. 7. Time evolution of (a) purity and (b) population with F = 1.6 x 10^^, for the LS 
case, atv^O (thick line),(- • •), ±3.9 ( ), ±4( ), ±5(- . - . -) and ±10 (thin line). 

Fig. 8. Upper panel: prr vs. the rescaled time. Lower panel: Purity vs. time. Here the LS 

case with v = ±5 is presented at F = (thick line), 1.6 x 10"*^ ( ), 1.6 x 10"^ (• • •), 

and 1.6 X 10-2 ( ). 

Fig. 9. Prr vs. the rescaled time in (a) and (c). Purity vs. the rescaled time in (b) and 
(d). The SDS case with F = 1.6 x 10"^ is presented. For (a) and (b), i; > 0: v^O (thick 

line), 2 ( ), 3 (■■■), 5(- .-.-), and 10 (thin line). For (c) and (d), v < 0: v^O (thick 

line), -2 ( ), -3 {■■■), -5(- . - . -), and -10 (thin line). Note the difference in the scale 

of ordinate in panels (a) and (c). 

Fig. 10. Prr vs. the rescaled time in (a) and (c). Purity vs. the rescaled time in (b) and 
(d). The WDS case with F = 1.6 x 10"^ is presented. For (a) and (b), v > 0; v = (thick 

line), 2 ( ), 3 (•••), 5(- . - . -), and 10 (thin line). For (c) and (d), t; < 0; ^; = (thick 

line), -2 ( ), -3 (•••), -5(- . - . -), and -10 (thin line). 

Fig. 11. Prr vs. the rescaled time in (a) and (c). Purity vs. the rescaled time in (b) and 
(d). The WDS case with v — 5 is chosen for (a) and (b) and the WDS case with v — —5 is 

chosen for (c) and (d). F = (thick line), 1.6 x 10"^ ( ), 1.6 x lO''^ (■■■), and 1.6 x 10"^ 

(— • — ). In panels (a) and (b) the F = and F = 1.6 x 10""^ curves overlap to the extent 
that they cannot be distinguished. 

Fig. 12. Prr vs. the rescaled time in (a) and (c). Purity vs. the rescaled time in (b) and 
(d). The SDS case with v = 5 is chosen for (a) and (b) and the SDS case with v — —5 is 
chosen for (c) and (d). F = (thick line), 1.6 x 10"^ ( ), 1.6 x 10"^ (•••), and 1.6 x 10"^ 
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Note the difference in tfie scale of ordinate in panels (a) and (c). 
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